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Abstract 

The role of viscous stress in heating of the fuel mixture in deflagration-to-detonation transi- 
tion in tubes is studied both analytically and numerically. The analytical theory is developed in 
the limit of low Mach number; it determines temperature distribution ahead of an accelerating 
flame with maximum achieved at the walls. The heating effects of viscous stress and the com- 
pression wave become comparable at sufhciently high values of the Mach number. In the case 
of relatively large Mach number, viscous heating is investigated by direct numerical simulations. 
The simulations were performed on the basis of compressible Navier-Stokes gas-dynamic equa- 
tions taking into account chemical kinetics. In agreement with the theory, viscous stress makes 
heating and explosion of the fuel mixture preferential at the walls. The explosion develops in an 
essentially multi-dimensional way, with fast spontaneous reaction spreading along the walls and 
pushing inclined shocks. Eventually, the combination of explosive reaction and shocks evolves into 
detonation. 
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I. INTRODUCTION 



Deflagration-to-detonation transition (DDT) is one of tfie most important and least un- 
derstood problems in combustion science. It is well known that a deflagration wave propa- 
gating from a closed tube end may spontaneously accelerate and trigger an explosion in the 
fresh fuel mixture, which goes over to detonation . The mechanism of DDT 

was described qualitatively already by Shelkin j^. According to the Shelkin explanation, 
thermal expansion in the burning process produces a strong flow in the unburnt fuel mixture, 
which becomes non-uniform because of non-slip at the tube walls. The non-uniform flow 
makes the flame curved, which increases the flame surface and hence the volumetric burning 
rate; it makes the flow even stronger and leads to flame acceleration. An accelerating flame 
pushes compression waves and shocks, the temperature in the fuel mixture increases, the 
reaction goes even faster, until explosion happens ahead of the flame front and evolves into 
detonation. 

For a long time it was a general belief that DDT is impossible without artificial turbu- 
lence generation. Since turbulence and turbulent burning are not so well-understood yet, 
it presents a major obstacle in developing a theory of flame acceleration and DDT. Only 
recently a constructive idea was suggested that DDT may be achieved even in a laminar flow 
in tubes with adiabatic walls P, |8|. Though it is difficult to obtain adiabatic walls in a real 
experiment, these conditions may be easily imitated in direct numerical simulations [3, Q]. 
Starting with this idea, the analytical theory of laminar flame acceleration was developed in 
[9], [lo| ; the theory was quantifled by extensive numerical simulations. The idea of laminar 
flame acceleration allowed also direct numerical simulations of DDT, e.g. see ,|12| and 
a recent review j^. The simulations demonstrated many interesting details of the DDT; at 
present it is not quite clear, which of them are intrinsic to DDT, making them the backbone 
of the process, and which are just supplementary effects of minor importance. One of these 
details is heating of the fuel mixture due to viscous stress, which makes the explosion pref- 
erential close to the walls. This is different from the classical Shelkin scenario of DDT with 
heating produced mainly by the compression/shock waves. 

A quantitative analytical theory of explosion triggering by an isentropic compression 



wave ahead of an accelerating flame was developed recently in 



13|. This theory disregarded 



influence of the viscous stress. However, an important role for viscous stress is in line (at 
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least, ideologically) with the model of the DDT, based on hydraulic resistance, see jl^ and 
references therein. This model had another shortcoming; it was one-dimensional and did 
not take into account Shelkin acceleration, which is an intrinsic part of DDT. Thus we came 
to the question, how strong the role of viscous stress is in heating of the fuel mixture and 
in DDT, taking into account the multi-dimensional Shelkin mechanism of this process. The 
purpose of the present paper is to answer this question. The paper considered only the 
geometry of smooth tubes, with no obstacles at the walls. 

We studied the role of viscous stress both analytically and numerically. The analytical 
study is based on the theory of laminar flame acceleration {q]. The theory works excellently 
for low Mach numbers, typical for the beginning of flame acceleration, about Ma = 10^^ — 
10~^. We found temperature profile in the flow ahead of the flame front with maximal 
temperature increase achieved at the walls. We show that, at low values of the Mach 
number, the role of viscous heating is extremely small, it scales as Ma^. By comparison, 
heating due to the compression wave scales as Ma. However, temperature at the walls grows 
faster in time than at the axis. According to the theoretical estimates, the heating effects of 
viscous stress and of the compression wave become comparable as soon as the Mach number 
(defined using total burning rate) reaches values above 0.6. Another interesting point is 
that heating due to viscous stress does not depend on the Reynolds number, at least as long 
as the Mach number is sufficiently low. When the Mach number becomes about 0.05 and 
higher, the theory works only qualitatively, providing order-of-magnitude estimates. 

In that range of Mach numbers we investigated the role of viscous heating by direct 
numerical simulations. The simulations were performed on the basis of compressible Navier- 
Stokes gas-dynamic equations taking into account chemical kinetics. It was again established 
that the fuel mixture is heated stronger at the walls than along the tube axis. In agreement 
with the analytical theory, we found that heating effects of viscous stress and compression 
wave become comparable at large values of the Mach number; viscous heating makes explo- 
sion of the fuel mixture preferential at the walls. The explosion develops in an essentially 
multi-dimensional way, with fast spontaneous reaction spreading along the walls and push- 
ing rather strong shocks inclined with respect to the tube axis. Eventually, the combination 
of explosive reaction and shocks evolves into detonation. 

The present paper is organized as follows: in Section II we develop the analytical theory 
of heating due to viscous stress; in Section III we describe the details of the numerical 
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simulations; in Section IV we present and discuss the results obtained; the paper is concluded 
with a short summary. 



II. THEORY OF HEATING DUE TO VISCOUS STRESS AHEAD OF ACCELER- 
ATING FLAMES 

We consider the same geometry as in jol : a laminar flame propagates in a two-dimensional 
(2-D) tube of half-width R with adiabatic walls and non-slip at the walls as shown schemat- 
ically in Fig. [H Burning matter expands as it passes the flame front; density ratio of the 
fuel mixture pf and the burnt gas ph is typically rather large, 9 = pj/pt = 5 — 8. Because 
of the thermal expansion, a flame propagating from the closed tube end pushes the unburnt 
fuel mixture and generates a flow. It was shown in jo] that the stream ahead of the flame 
may be closely approximated by a plane-parallel flow along the walls u = ezUz{x,t). 

The normal planar flame velocity Uf provides a natural scaling of velocity dimension for 
the problem. However, the total burning rate Uyu is different from Uf-. it shows how much 
fuel mixture is consumed per unit time by the whole flame front and how much energy is 
produced. In the standard 2-D model of an inflnitely thin flame front, the relative increase 
in the burning rate is simply equal to the increase in the total length of the flame front. The 
propagating flame generates a flow, which is not uniform because of non-slip at the walls. 
Friction retards the gas close to the walls, while flow velocity at the tube axis is larger than 
the average one. The non-uniform velocity proflle distorts the flame shape, which leads to 
increase in Uw and hence to flame acceleration. At sufficiently small values of the Mach 
number, neglecting the starting transitional period, the flame has been shown to accelerate 
exponentially [9| 

U^/Uf = exp (crUft/R) . (1) 



For 2-D flow, the dimensionless acceleration rate a was found to be 
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where the value Re = UfR/u plays the role of the Reynolds number for the problem. For 
large values of the Reynolds number, Eq. ([2]) predicts decrease of the acceleration rate 

a = Q^Re. (3) 
4 
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In [9j, the Navier-Stokes equation was solved as 



.^ , cosh u — cosh (ux/i?) 

VC/; = (9-l)exp(.C/,t/fl) J^^_J^\ (4) 



where ^ = \J oRe. For large values of the Reynolds number we have fi = Q, and the flow 
resembles qualitatively a combination of two boundary layers at the walls separated by the 
main almost uniform stream, see Fig. O Paper demonstrated a very good agreement of 
the theory with direct numerical simulations of acceleration flames. 

Here we find the temperature increase in the fuel mixture ahead of the accelerating flame 



using the results above. Similar to 



9j, we adopt the limit of incompressible flow, which 



holds in the case of a sufficiently small Mach number. In that case the equation of thermal 



conduction in the fuel mixture is 



15| 



^ + u.vr = ,v^r + — + . (5) 

In the case of a shear flow u = e^Uzix, t), with T = T{x, t), Eq. ([5]) reduces to 



dt ^dx^ ^ Cp[dx ) ■ 

We introduce standard scalings R, Uf and R/Uf as units of length, velocity and time. In 
that case we work with the dimensionless values {r];^) = {x;z)/R, t = Uft/R, w = u^/Uf, 
and the dimensionless temperature = T/Tq — I. Taking into account the relation for sound 
speed in a politropic gas c^q = (7 ~ l)Cp^o with heat capacity Cp and adiabatic exponent 
7, we rewrite Eq. dH]) as 

_ 1 d'^d Mai (dwV 

~ PrRe dr]^ ^ ~ ) ' 

where Moq = Uf/cgo is the Mach number deflned with the help of the planar flame velocity. 
The last term in Eq. (JTj) is specifled by the theory with 

(^)^(e-l)Vexp(2..) ^^^^^;^^'^t.,). - (S) 

We evaluate the possible effect of heating keeping in mind the limit of large Reynolds number 
with cr = Q'^/Re, /i = 0, see Eq. ([3]). Maximal heating happens close to the wall. Due to 
symmetry, we may consider only the wall at rj = 1, where sinh(/i?7) ^ cosh(/i?7) ^ |exp(/i?7). 
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Such an approximation holds with exponential accuracy exp(— 0) << 1. In that case we 
have close to the wall (at 1 — r/ << 1) 

dw\ _ (6 - l)'^4g^p [2crr + 2(1(7] - 1)] , (9) 



drj J (/i — 1)2 
and dw/drj ^ at the channel axis rj = 0. Then Eq. ([7]) reduces to 

1 9¥ , Ma2(e-p2 



A particular solution to Eq. (JTOi) may be written as oc exp [2crr + 2/1(77 ~ 1)]) which leads 

to ^ 

^ _ (7-l)Pr(9-l) /x^ ^^, _ ^ 

2 (2 - Pr) (/I - 1)' ° ^ A-v/ ;j v ; 



Here we take into account that /i = yoRe. A complete solution to Eq. (ITOl) is a superposition 
of the particular solution Eq. (ITTl) and the solution to an auxiliary equation 

dd 1 (92^9 



9r PrRe dr]"^ 
which also grows in time as oc exp(2(Tr) so that 



(12) 



(9¥ 

2aPrRe^ = — . (13) 
Solution to the auxiliary Eq. f|T3l) may be found as 

i^a = [Ai exp{kr]) + ^2 exp(-A;?7)] exp(2crr), (14) 
with unknown amplitudes Ai, A2 and 



k = V2aPrRe = V2Prfi. (15) 



The complete solution to Eq. ( fTOl) is a superposition of and -^a, Eq- ( pAl) and Eq. (|T4l) . 
satisfying the adiabatic boundary condition at the wall and the symmetry condition at the 
channel axis 

(9^9 

- = 0, v = 0,l. (16) 

Since "i^p is exponentially small at the axis r/ = 0, we find Ai = A^. Within the same 
exponential accuracy, the decaying term exp(— /cr/) is negligible close to the wall. Thus 
we find at the wall rj = 1 

kA, exp = (7 - 1) ^2 _ ^ _ [^^ Mai (17) 



which determines the complete solution close to the wall as 



2 ,,2 



(7- l)Pr(0- 1)V 



■Mai exp(2crr) 



2fi 
T 



2 (2 - Pr) in - 1 
The temperature increase at the walls rj = 1 is 

(7-l)(0-l)V^ ViV 



exp [k{ri — 1)] — exp [2fi{ri — 1)] 



(18) 



-.Malexp{2aT). 



(19) 



2 (id -if V2 + VP^^ 
Finally, as an evaluation, we consider the limit of large Reynolds number a = Q^/Re, fi = Q, 



which leads to the temperature increase 



Ma5exp(2crr) 



--U'Jcl. (20) 



2 ~ V2 + VP^ 2 V2 + VPr 

It is interesting that viscous heating at the wall, Eq. fl20l) . does not depend on the Reynolds 
number. Combination Ma^ = Maoexp((Tr) = Uw/csq determines the running value of the 
Mach number for an accelerating flame, which is different from the initial value Moq = 
Uf/csQ. As an example, we take 6 = 8, Pr = 0.7, 7 = 1.4, and find the temperature 
increase due to viscous stress {} ~ 4.8Ma^ = 4.8[/^/c^q. 

We compare this result with the heating due to the compression wave found in 

1 + 7? 



1 + ^^^ — ^(6- l)Ma, 



131 ] as 
(21) 



Note that the designation d has a different meaning in the present calculations and in 
In the case of small Mach numbers we find 
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t9 = (7 - i)(e - i)Ma^ = (7 - i)(e - 1) 



(22) 



Using the same numbers as before we obtain a temperature increase in a compression wave 
d = 2.8Uw/cso- Thus, for = 8, Pr = 0.7, 7 = 1.4, we should expect a comparable 
temperature increase due to the compression wave and due to viscous stress at the Mach 
number about Ma^ = U^/cso ~ 0.6, see Eqs. ( I20l) . (!22l) . Besides, heating due to the 
compression wave is mostly a ID process equally strong at the walls and at the channel axis. 
As both effects work together, we find preferential heating at the walls. However, we have 
to remember that the theory of flame acceleration 9| works quantitatively well only for a 
noticeably subsonic flow. Ma = 10^^ — 10^^. At such small Mach numbers heating due to 
viscous stress is negligible, since it scales as ?9 oc Ma^, while heating in a compression wave 



has another scahng Ma oc Ma^. At Mach numbers above 0.05 the theory provides only 
order-of-magnitude estimates. Direct numerical simulations of the present paper show that, 
at Mach numbers above 0.05, flame acceleration goes much slower than in the exponential 
regime Eq. ([T]). The critical value Ma^ = Uy^/cgo ~ 0.6 obtained above is also an order-of- 
magnitude estimate; it is beyond the rigorous limits of the developed theory. Still, the result 
Eq. (120|) may be useful even at relatively high values of the Mach number. For comparison, 
temperature increase due to viscous stress in the stationary boundary layer produced by a 



uniform flow Ua is 

AT 7-1 , — f/2 ^ ^ 

Here and c^a are the temperature and the sound speed in the uniform flow far away 

from the wall. In the present problem flow velocity in the main stream (e.g. at the axis) 

plays qualitatively the same role as Va- Rewriting Eq. ( l20l) in the same form, we flnd 

AT 7-1 /P^ Ul , , 
= 1= — . (24) 

Here we have taken into account that velocity at the axis ahead of an accelerating flame is 
related to the total burning rate as Ua = QUw, see Eq. @. The difference between Eqs. f l23|) 
and (12^ is then simply in the numerical factor \/2 + y/Pr, which is about 2.25 for Pr = 0.7. 
An exponentially accelerating flame produces about half the temperature increase at the 
wall compared to that from the stationary boundary layer. Thus one would expect that the 
temperature increase at Ma^ = 0.05 and above is somewhere between predictions of Eqs. 
f l23|l and flMl) . In order to investigate the role of viscous heating at relatively large values 
of the Mach number we performed direct numerical simulations. We stress once more that 
both formulas ( l23l) . ( !24l) predict viscous heating independent of the Reynolds number. This 
conclusion is of great help in numerical simulations, since it allows investigating the process 
even for rather low values of the tube width. 

III. BASIC EQUATIONS OF THE DIRECT NUMERICAL SIMULATIONS 

We performed direct numerical simulations of the 2-D hydrodynamic and combustion 
equations including transport processes (thermal conduction, diffusion, viscosity) and chem- 
ical kinetics with an Arrhenius reaction. The equations read 
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(26) 



dt V 2^^'^' ) dx'- \P^^ ^ -pUiUiUj + qj - u^ij J = 0, (27) 

I (pn + I- {pn.Y - Ag:) . e.p i-E/R,T) . (28) 

where Y is the mass fraction of the fuel mixture, e = QY + CyT is the internal energy, 
h = QY + CpT is the enthalpy, Q is the energy release in the reaction, Cy, Cp are the heat 
capacities at constant volume and pressure. It was assumed that the heat capacities do not 
depend on the chemical reaction. The stress tensor and the energy diffusion vector qj 
take the form 

/ dui dui 2 duk r \ 

where /i = pz/ is the dynamic viscosity, Pr and Sc are the Prandtl and Schmidt numbers, 
respectively. To avoid the Zeldovich (thermal-diffusion) instability we take the unit Lewis 
number Le = Pr/Sc = 1, with Pr = Sc = 0.75. The dynamical viscosity is /i = 1.7 x 
10~^ Ns/m'^. The fuel-air mixture and burnt matter are assumed to be a perfect gas with 
a constant molar mass m = 2.9 x 10~'^kg/mol, with Cy = 5-Rp/2m, Cp = 7Rp/2'm, where 
Rp ~ 8.31J/ (mo/ ■ K) is the perfect gas constant. The adiabatic index is 7 = Cp/Cy = 1-4. 
The equation of state is 

P = pRpT/m. (31) 

We used the initial temperature of the fuel mixture Tj = 300K, initial pressure Pj = lO^Pa, 
and initial expansion ratio Gq = 8. The initial Mach number is Moq = Uj/cgQ = 0.04. Eq. 
(I28l) describes a single irreversible reaction of the first order, where the reaction rate obeys 
the Arrhenius law with the activation energy Ea and the factor of time dimension th- In the 
simulations, the value of Uj is determined by the choice of thermal-chemical parameters such 
as Ea, Tji, and the energy release in the reaction Q = (Bq — l)CpTf. In order to have better 
resolution of the reaction zone, we took Ea/RpTj = 32 in the simulations. Note that the 



process of viscous heating addressed in our paper does not depend on the activation energy. 
The factor was adjusted to obtain a particular value of the planar flame velocity f// by 
solving the associated eigen-value problem. The flame thickness is conventionally defined as 

where pj = l.lQkg/m? is the initial mixture density. However, we would like to stress that 
the value (1321) is just a thermal-chemical parameter of length dimension in the problem; the 
characteristic size of the burning zone rnay be an order of magnitude larger [17l]. 

In the studies of flame acceleration [oljIlO], the main parameter of simulations was the tube 
width 2i?, which controlled the Reynolds number. In the present paper we are interested 
in thermal heating due to viscous stress. According to the present theory and the classical 
results [l6|, viscous heating does not depend on the Reynolds number (or this dependence is 
extremely weak). For this reason, to study the effect, it is sufficient to perform simulations 
only for one channel width. The channel width was set to 2R = 40L/, which implied the 
Reynolds number related to the planar flame velocity Re = UfR/v = R/LfPr = 26.7. The 
Reynolds number related to the flow Re = {uz)'2R/v may be larger by several orders of 
magnitude due to flame acceleration and thermal expansion of the burning gas. We took 
non-slip and adiabatic boundary conditions at the channel walls: 

u = 0, n ■ VT = 0. (33) 

where n is a normal vector at the wall. At the open channel end non-reflecting boundary 
conditions are used; the boundary conditions were tested in [l^]. We used two types of 
initial conditions: 1) a hemispherical flame "ignited" at the channel axis at the closed end 
of the channel, or 2) a planar flame front. 

The internal structure of the flame imitated the analytical solution by Zel'dovich and 
Frank-Kamenetskii [31]. Particularly, for the hemispherical case it was given by 

T = + Tj(e - 1) exp (^-Va;2 + z'^/L^^ , if z^ + > r). (34) 

T = GTj, if + < r). (35) 

F=^^^, P = P,, u = 0. (36) 
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Here Vf is the radius of initial flame ball at the closed end of the channel. For an initially 
planar flame shape Eq. flMl) is replaced by 



T = Tf + Tf{e -l)exp[- 



(37) 



if 2; > Zo, with Eq. fl35|) for z < zq [zq denotes the initial flame position). 

We simulated flame dynamics and DDT using a 2-D Cartesian Navier-Stokes solver with 



chemical reactions developed at Volvo Aero. T 



le code is based on the cell-centered flnite- 



volume scheme developed by L.-E. Eriksson 18 



191]. The numerical method has proved to 



be robust for modelling different kinds of complex reacting flows. The TVD-limiting is used 



where applicable to prevent overshoots in flow fleld properties in regions o 
The code has been validated by solving various hydrodynamic problems 



high gradients. 



18|, lia, 120[, and 



was utilized successfully in simu 



0,0, 



12 



17 



21 



22 



23 



24 



25 



ations of laminar flames at different conditions of burning 



26|. 



The grids used in present simulations were non-uniform with higher resolution around 
the flame front and pressure waves and lower resolution in other parts of the calculation 
domain. Grid resolution varied gradually from one zone to another along the tube. Using 
this method we avoided processing a too large number of cells. Square cells of the size 0.2L/ 
were used in the flame domain to ensure isotropic propagation of the curved flame in x and 
y directions. In |l7| we demonstrated that such a grid is sufficient to resolve even a strongly 



curved flame front. As an illustration, in Fig. 3(a) we present the proflles of the main 



variables inside the flame front: the temperature distribution (T — Tf)/(Th — Tf), the mass 
fraction Y and the reaction rate A = [pY/rfi) exp {—Ea/RpT) scaled by its maximal value. 
We observe that the thermal-chemical parameter Lf deflned by Eq. ( |32l) is much smaller 



than the real flame thickness. According to Fig. 3(a), the characteristic length scale of the 



temperature proflle is about (4 — 5)Lf, while the width of the active reaction zone is about 
Lf. Taking the grid size of 0.2Lf in ^-direction we obtained about flve grid points inside 
the active reaction zone and 20-25 grid points inside the flame. 



A sketch of the calculation grid is shown in Fig. 3(b) The calculation grid is modifled 
periodically in the course of calculation, so that the sizes and positions of all domains in Fig. 



3(b) change dynamically, following the flame front and pressure waves. Thus the leading 



compression wave can never reach the open tube end during a calculation. Cell size was 
0.4L/ around the front of the leading pressure wave. Splines of the third order are used 
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for re-interpolation of the flow variables during periodic grid reconstruction to preserve the 
second order accuracy of the numerical scheme. We performed several test simulation runs 
with different resolution between 0.1 and 0.2Lf in the flame domain. The tests demonstrated 
that the difference in the flame velocity, pre-detonation time and distance does not exceed 
10% for different resolutions. 



IV. SIMULATION RESULTS 

We performed direct numerical simulation of flame dynamics and detonation triggering 
in a tube of width 2R = 40Lj with non-slip walls. The flame propagated from the closed 
tube end to the open one; a semi-infinite tube was imitated by using non-reflecting boundary 
conditions at the open end and by adjusting the mesh to the leading pressure wave. We 
started the simulations either with a hemispherical flame ignited at the tube axis or with a 
planar flame front. Fig. H] shows velocity of the combustion front in the laboratory reference 
frame of the tube, Uiab/Uf, in the case of hemispherical ignition; the velocity was calculated 
as the average over the reaction front. Flame dynamics involves several velocity parameters, 
which have quite different physical meaning and may differ quantitatively by one or two 
orders of magnitude. First, we have a planar flame velocity Uf, which characterizes local 
normal velocity of the flame front relative to the moving fuel mixture (at least in the model 
of an infinitely thin front). The value Uf works as a standard velocity scaling in the problems 
of flame dynamics. To be accurate, Uf depends on thermal and chemical parameters of the 
fuel mixture, which vary in the process of detonation triggering. Here Uf designates the 
initial planar flame velocity. Second, we have the burning rate U^, which specifies total 
amount of fuel mixture burnt per unit time (see also Section II). The burning rate f/^ may 
exceed Uf hj several orders of magnitude in the process of flame acceleration. Still, the 
value Uw does not describe flame propagation in any particular reference frame except for 
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26 



27 



the case of stationary or statistically-stationary curved flames like in 

Flame propagating from the closed end pushes a strong flow in the fuel mixture. In the 
incompressible case, average velocity of the flow is (m^) = (0 — l)f/tt,, where averaging is 
performed over the channel cross-section. The flow drifts the flame, and the drift velocity 
is different for different parts of the flame front. For this reason, the burning rate U^ is 
also different from flame velocity in the laboratory reference frame, with respect to the tube 
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walls, Uiab- We stress that Fig. H] presents flame velocity in the laboratory reference frame, 
Uiab, not the burning rate 1/^. In the case of an essentially subsonic flame, we have a simple 
relation Uiab = ^U^j, which may also provide an estimate in the case of flnite Mach numbers. 
Different velocity parameters imply different Mach numbers as well, which we designate by 
Moq = Uf/cso, Ma^ = Uu,/csa, Maiab = Uiab/cga', label "a" refers to the position at the 
channel axis just ahead of the flame front. Taking a rather low initial Mach number related 
to the planar flame velocity Mao = 0.04 with O = 8, we have immediately Maiab = 0.32, 
which is not so low. Soon after ignition, in a short time about (0.2 — 0.3)R/Uj, burning 
rate increases by order of magnitude in the process of precursor acceleration 29|, l30|]. As a 
result, quite quickly the Mach number Maiab becomes comparable to unity. We also point 
out that Mqwi Maiab take into account variations of the sound speed in the fuel mixture in 



time and space; scaled sound speed Csa/Uf is shown by the dashed line in Fig. 4(a^ 



Transition from one combustion regime to another (flame-explosion-detonation) may be 
clearly seen by drastic changes in the front velocity. In the flame regime, we can see front 
acceleration, which involves two different physical mechanisms. First, we observe a short 
period of precursor acceleration, related to the hemispherical geometry of flame ignition. 
This type of acceleration is possible even in the case of friction-free walls. Precursor ac- 
celeration goes noticeably faster than the Shelkin mechanism, but this process is limited in 
time. Precursor acceleration for a cylindrical tube was investigated experimentally in 29 1 
and numerically in |30|]. The paper 30 1 presented also the analytical theory of the process. 
In the present geometry of a 2-D channel the precursor acceleration is somewhat different 
from the cylindrical case; still, the basic features of the process remain the same. 

The precursor acceleration is followed by acceleration due to the Shelkin mechanism 
produced by thermal expansion and non-slip at the walls. The theory of the Shelkin ac- 
celeration in incompressible flows has been developed in [o], lo|. However, in the present 
case Mach number Maiab is comparable to unity or above, and the incompressible theory 



of flame acceleration describes the process only qualitatively. A dashed line in Fig. 4(a^ 



shows the exponential growth predicted by the incompressible theory Eq. ([T]); it goes well 
above the numerical results as soon as the precursor acceleration is over. Investigation of 
the compressibility effects on flame acceleration is beyond the scope of the present paper; 



this is work for the future. Still, as a preliminary result. Fig. 4(b) shows the flame veloc- 
ity in log-log variables. In order to elucidate the transition effects related to a particular 
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geometry of flame ignition, we have performed two simulation runs: that of a flame ignited 
at the tube axis and one for a flame ignited as a planar front. In the flrst case we have 
precursor acceleration of the flame front, while in the second case this effect is missing. As 
we can see, both simulation runs tend to the same asymptotic regime of flame acceleration 
with Uiab oc (Uft/R)"', where n ^ 0.36. More detailed investigation of the process will be 
presented elsewhere. It is well known that flame is subsonic; still, this rule applies to the 



burning rate Uw, not to Uiab- In Fig. 4(b), Uiab crosses the sonic line at the time instant 
about Uft/R = 30. This time instant is often called run-up time to supersonic flames |31 |. 

An accelerating flame produces a flow in the fuel mixture, which is non-uniform because 
of the non-slip boundary conditions. A characteristic velocity proflle is shown in Fig. O 
the results of numerical simulations are compared to the predictions of the incompressible 
theory [q]. The velocity proflle is calculated at different time instants just ahead of the flame 
front and scaled by the maximum velocity value at the tube axis. Both the theory and the 
present simulations demonstrate the same qualitative structure of the flow, which consists 
of an almost uniform stream close to the axis and transitional (boundary) layers close to 
the walls. The velocity proflle is in a good agreement with the incompressible case, not 
only qualitative, but even quantitative. A non-uniform viscous flow implies stress, which is 
stronger close to the channel walls. As a consequence, one should expect additional heating 
of the fuel mixture due to the viscous stress, as discussed in Section II. In this Section, our 
main goal is to investigate the process of viscous heating numerically and to compare it to 
the theory of Section II. 

The temperature increase in the fuel mixture just ahead of the flame front is shown in Fig. 
|6] as a function of time and at two positions: at the tube axis and at the channel wall. We 
can see that, within a short time after ignition Uft/R = 0.2 — 0.3, the temperature increase 
is already different from zero, being about AT/Tq ^ 0.25. This happens because of the 
leading shock wave, pushed by the flame in the process of precursor acceleration. The flow 
velocity behind the shock is about Auz/Uj ~ — 1, and thermal expansion is typically quite 
strong. Therefore, in spite of a small initial Mach number MaQ = 0.04, the flow velocity is 
not so small in comparison with sound speed even for a planar flame front, Am^/cso = 0.28. 
Precursor acceleration makes this value even larger. At the beginning, temperature increase 
at the walls and at the tube axis is practically the same, which illustrates almost ID structure 
of the compression wave. As time goes on, the flame accelerates, the Mach number of the 
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flow increases and heating due to viscous stress becomes noticeable in agreement with the 
theory of Section II. By the time of explosion, Uft/R ^ 4.38, difference between temperature 
at the axis and at the walls is noticeable; still, this difference is not as strong as predicted 
by the theory. In order to compare the numerical simulations and the theory quantitatively, 
in Fig. [7] we plot temperature increase versus Mach number Maiab- In that case we separate 
heating in the compression wave (at the axis), and the effect of viscous stress as temperature 
difference between the wall and the axis, T^aii — Taxis- Heating in the compression wave at 
the axis is shown by squares. Temperature at the axis increases almost linearly with the 
Mach number, in agreement with the theory, Eq. ( 12T|) . presented by the solid line. The 
theoretical prediction holds for a ID isentropic flow; small deviations of the theory and 
simulations indicate that the flow is not completely ID and isentropic. Two other solid 
lines show theoretical predictions for temperature increase due to viscous stress ahead of an 
exponentially accelerating flame, Eq. (12^ . and in a stationary boundary layer, Eq. (!23ll . 
These two formulas present two limiting cases of strong acceleration and no acceleration at 
all. As we pointed out, at relatively high values of the Mach number, the flame acceleration 
is much slower than exponential. By this reason it is natural to expect that the numerical 
points should be somewhere in between the predictions of Eqs. fl23|) and fl24l) . Indeed, Fig. [7] 
shows that this is the case; the numerical results come rather close to that of Eq. (12^ . being 
a little above the theoretical curve. Only at the end of the simulation run, temperature at the 
wall increases drastically, which indicates the beginning of the explosion. It is also interesting 
to look at the temperature distribution in the channel. Fig. [8] compares temperature profiles 
at different distances ahead of the flame front at the time instant Uft/R = 3.43 to the 
theoretical result Eq. (fT8|) . The theoretical predictions look quite similar to the numerical 
results: all the time we observe the effect of viscous heating well-localized at the walls. Still, 
the theory predicts a self-similar temperature profile independent of position, since it was 
developed for an exponentially accelerating flame in an incompressible flow. In the numerical 
simulations flame acceleration is not exponential, the flow is not incompressible, and the 
temperature profile changes in space and time. As we can see in Fig. [8], the transitional zone 
of viscous heating is wider close to the flame front, and it becomes more localized far away 
from the flame. 

As we can see, heating due to viscous stress is not large in comparison with that in a 
compression wave, and, in principle, explosion theoretically could be possible even without 
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viscous heating. Fig. [9] presents temperature field snapshots in the area of the accelerating 
flame from beginning of the explosion till detonation. Prior to the explosion, we can see light 
blue strips along the walls indicating regions of viscous heating. Still, Arrhenius reactions are 
extremely sensitive to temperature; even this slight non-uniformity influences strongly the 
structure of the explosion. The explosion may be observed already at the second snapshot: 
it develops in the form of two narrow tongues spreading fast along the tube walls from the 
flame front to the fuel mixture. Temperature in the explosion zone is noticeably higher 
than it was behind the flame front. The tongues grow in length and width, until the original 
flame front becomes engulfed in the explosion. However, even at that time the explosion front 
remains essentially multi-dimensional: a large amount of the unburned matter separates two 
explosion tongues. The explosion pushes inclined pressure waves, well seen in Fig. M by the 
light blue color. Intersection of the pressure waves increases temperature at the channel axis, 
reduces the reaction time, until a bridge is formed between two tongues of the explosion. 
As two tongues meet, they soon develop into detonation, which can be seen in the last two 
snapshots of Fig. [9l Meeting of two explosion tongues may be also seen in Fig. |H as an 
extremely sharp peak in the velocity dependence. Strictly speaking, this peak is infinitely 
sharp and high, since it corresponds to a singularity of pocket formation. In Fig. |U we cut 
down the peak to make the plot more illustrative. In the present calculations, DDT distance 
is 330i?, but this value is, of course, sensitive to a particular chemical kinetics. 



V. SUMMARY 



The purpose of the present paper was to investigate the role of viscous heating in DDT. 
We solved this problem both analytically and by direct numerical simulations. The ana- 
lytical theory was developed within the approximation of incompressible flow, which holds 
at sufficiently low values of the Mach number. On the contrary, the numerical simulations 
were performed for relatively large Mach numbers. Both theory and numerical simulations 
demonstrated that the role of viscous heating increases with the Mach number. Viscous 
stress makes temperature higher at the channel walls in comparison with the axis. Temper- 
ature profile in the channel is given by Eq. ( fTSl) . temperature increase at the wall due to 
viscous stress is determined by Eq. fl2^ . The theoretical results agree rather well with the 
numerical simulations, though, when making a comparison, one has to take into account 
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much slower flame acceleration at larger values of the Mach number. Even just before the 
explosion, additional temperature increase due to viscous heating is noticeably smaller than 
the temperature increase in the compression wave. Thus, viscous heating should be treated 
only as a supplementary process in DDT. In principle, DDT could be possible with heating 
provided by the compression wave only in line with the Shelkin scenario Still, 
neglecting viscous heating in a theoretical model like [3] may lead to considerable disagree- 
ment in the time and distance of detonation triggering predicted by the theory and obtained 
in experiments/numerical simulations. Viscous heating makes explosion triggering easier 
and faster because of extreme sensitivity of Arrhenius reactions to temperature variations. 
In a certain sense, it plays the same role as hot spots: it is not necessary, but helpful. Non- 
uniform temperature distribution leads to an essentially multi-dimensional structure of the 
explosion, which pushes inclined pressure waves and eventually develops into detonation. 
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FIG. 1: Flame propagating in a tube with non-slip walls. 
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FIG. 2: The velocity profile Uz scaled by the amplitude Umax- The solid line shows the theoretical 
result Eq. The markers correspond to the simulation results of [9] for Re = 50, Pr = 0.5 

at the distances lOR; 20R; 35R (circles, triangles and crosses) from the flame at the time instance 
Uft/R = 1.15. The arrows illustrate the direction of the flow. 
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FIG. 3: Illustration of the grid, (a) Profiles of the scaled temperature (T — Tf) / {Th — Tf), the mass 
fraction Y and the reaction rate A = (pY/tr) exp (—Ea/RpT) scaled by its maximal value within 
the flame front, (b) The adaptive grid used in numerical simulations. 



21 








\ 

00 jfTSfiw^er mf:^d^^m 



0.01 



0.1 



10 



(b) 

FIG. 4: (a) Scaled velocity of the reaction front (fiame-explosion-detonation) in the laboratory 
reference frame, Uiah/Uj. One dashed line presents theoretical predictions for the incompressible 
flow. The other dashed line shows scaled sound speed, Cga/Uf. (b) The same plot in log-log 
variables for a flame ignited at the channel axis (with precursor acceleration) and as a planar front 
(no precursor acceleration). 
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FIG. 5: The velocity profile Uz scaled by the amplitude Umax- The solid line shows the theoretical 
result Eq. dH). The markers correspond to the present simulation results just ahead of the flame 
front at the time instants Uft/R = 0.92; 1.74; 3.43. 
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FIG. 6: Temperature increase at the wall and at the channel axis versus time. 
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FIG. 7: Temperature increase versus Mach number in tlie laboratory reference frame. Solid lines 
show theoretical predictions for temperature increase in the compression wave at the axis, Eq. 
(pij) . and additional increase due to viscous stress at the wall, Eqs. ([25]) . ([2^ . The dashed lines 
and the markers show respective numerical results. 
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FIG. 8: Scaled temperature profile (T — Taxis) /{Twaii — Taxis)- The solid line shows the theoret- 
ical result Eq. (jlSp . The markers correspond to the present simulation results at the distances 
(6; 425; 930)Lj ahead of the flame front at the time instant Uft/R = 3.43. 
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FIG. 9: Temperature field during DDT for the time period from Uft/R = 4.38 to Uft/R = 4.58 
witli equal time intervals between the snapshots. 
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